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ABSTRACT 


The  laminar  wake  of  a  re-entry  vehicle  is  simulated  by  a  system 
of  cylindrical  plasma  shells.  l'he  radar  cross-section  per  unit  length  is 
found,  taking  into  account  the  energy  losses  due  to  collisions.  Numerical 
results  await  completion  of  a  computer  program. 
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INTRODUCTION 


It  has  been  observed  that  a  radar  echo  is  received  from  the  wake  of 
a  vehicle  upon  re-entry  into  the  atmosphere,  and  that  tins  echo  represents 
a  considerable  enhancement  of  the  total  radar  return  from  the  re-entry 
event.  The  purpose  of  this  paper  is  to  provide  a  theory  which,  under 
certain  conditions,  will  predict  the  return  from  such  a  wake. 

The  laminar  wake  is  assumed  to  be  of  such  form  that  it  can  be 
simulated  by  cylindrical  shells  concentric  about  the  axis  of  the  vehicle, 
each  shell  being  a  homogeneous,  isotropic  medium.  The  media  are 
assumed  to  be  weakly  ionized  plasmas  containing  no  true  charge;  i.  e.  , 
the  net  charge  in  any  finite  volume  is  zero.  A  complex  permittivity  is 
expressed  in  terms  of  plasma,  collision  and  excitation  frequencies  in 
order  that  a  symmetrical  wave  equation  may  be  used  to  describe  the 
electromagnetic  fields.  Initially,  a  plane  wave  of  arbitrary  orientation 
with  respect  to  the  wake,  is  assumed  to  be  incident  on  the  outer  shell. 
Complexity  of  the  equations  requires  that  only  normal  incidence  be 
considered  at  this  writing.  Solutions  for  the  fields  in  the  shells  are 
found,  and  the  resulting  scattered  wave  is  found  by  matching  the  fields  at 
the  various  boundaries.  Once  the  scattered  field  is  determined  at  the 
surface  of  the  outer  shell,  the  field  at  an  arbitrary  observation  point 
may  be  found  by  using  the  Ki rchhoff-lluygetis  principle  for  vector  waves 
as  derived  by  Stratton  (5).  The  radar  cross-section  is  then  calculated. 
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ANALYSIS 


Thu  fields  in  the  various  regions  must  satisfy  Maxwell's  equations. 
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The  constitutive  equations  are: 

D-  eE  (5) 

B  • /I  U  (6) 

J  -  0-£  (7) 


Since  we  have  assumed  constant  permittivity,  charge-free  regions, 


v  •  £  =*  o 


(8) 


Assuming  an  exp*  ()(»)t)  time  dependence  for  the  field  vector,  (2)  may  be 
written 
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v lu* rt*  cjp  is  tiu*  plasm.*  frequency*  y'  is  the  collision  frequency  and  oJ  ,  tin* 
*  i  tat  toil  Irequent  y  (S).  I  lien 
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1-  l  oin  (I),  (h) ,  and  (H*)*  we  w  rite 
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latlin^  £  yU  ~  K  I  rum  (K),  we  may  write 
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Usiiu!  the  s.iim*  method*  it  may  be  shown  lh.it 


V*  H  +  K 2  H  -  C 
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With  the  assumed  time  dependence,  the  square  rout  of  must  be  extracted 

such  that  the  sign  of  the  imaginary  part  is  negative.  The  Laplacian  of  a 
vector  is  defined  to  be  the  sum  of  the  Laplacian's  of  the  rectangular  com¬ 
ponents  of  the  vector.  Therefore,  either  (11)  or  (12)  may  be  written 

4t 

where  CL  =  x,  y,  z  and  is  a  unit  vector.  Due  to  the  orthogonal  proper¬ 
ties  of  the  components,  each  scalar  component  must  satisfy  the  equation 
independently  so  that 

+  <u> 

In  a  transformation  of  co-ordinates  from  rectangular  to  cylindrical, 
the  z  component  is  unaffected.  Therefore,  a  solution  of  (13)  is  a  solution 
for  the  z  component  of  the  E  or  H  field  in  cylindrical  co-ordinates.  In 
cylindrical  co-ordinates  (1  3)  becomes 

»!£  + 1 Ut  ♦  •  »1±  ♦  ill  *  kU  =  o 

9r*  y  Sr  r  SB 1  9**  (14) 


The  variables  may  be  separated  so  that 

<l>  -  0(0)  R(r)T  (*) 

is  a  solution  of  (14).  The  equation  is  then  satisfied  for 
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where  /.  (r/K2  -  h^)  is  a  solution  of  the  Bessel  equation  and  h  is  an  arbitrary 

constant. 

The  following  are  particular  solutions  of  the  Bessel  equation,  and 
their  important  characteristics  (Z).  The  symbol ^  will  represent  \j -  h" 

I. 
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If  n  is  replaced  by  -n  in  (16),  the  equation  is  unaltered  and  .1  (jlr)  is  also 
a  solution.  However*  fur  integral  n,  the  two  solutions  are  not  litiearlx 
independent.  Jn  </r)  and  J_n  (^r)  .ire  known  as  Bessel  functions  of  the 
first  kind.  They  are  the  only  solutions  o!  (lb),  whirl  remain  finite  for 
r  b  fl, 

II. 

W-  dr  nir  [j»  u )  n  7r  J  n  ( $  r  )  ] 


This  is  known  as  a  Neumann  function  or  Bessel  function  id  the  second  kind. 
For  integral  n,  this  function  becomes  indeterminate  but  max  l»e  evaluated  h. 
L'lfospital's  rule. 


The  Bessel  functions  of  the  first  and  second  kind  converge  for  ail 
values  of  the  argument,  real  or  complex,  except  as  noted  at  the  origin. 
However,  for  large  arguments  the  convergence  is  slow  and  asymptotic 
forms  must  be  used. 

III. 

h/Vaj-  j,ca)+  (*)  09, 

(fir)  m  Jfy  (At)  -  j  fi*  (fir)  (20> 

These  functions  are  known  as  Hankei  functions  of  the  first  and  second 
kind  respectively.  These  functions  vanish  at  infinity.  With  the  assumed 
time  dependence  exp  (j*>t),  and  Hn^  (/&)  describe  waves 

traveling  radially  inward  and  outward  respectively. 

Consider  a  set  of  three  concentric  cylinders,  the  outer  cylinder 
bounded  by  infinity  (Figure  1).  The  wave  functions  for  the  three  regions 
may  now  be  determined.  Region  I  contains  the  origin  and  thus  /.  (^r)  must 

be  a  Bessel  function  of  the  first  kind. 

*X  -  ait  *  (/*r) (2|) 

Region  II  will  contain  waves  transmitted  from  Region  III  and  waves  reflected 
from  the  surface  at  Region  I,  Therefore,  a  combination  of  Hankei  functions 
of  the  first  and  second  kind  is  appropriate. 

[  lv  H„W  (A)  *  C,  <22> 


7 


Region  III  contains  .«n  incident  wave  yet  to  be  specified  and  a  scattered  wave 
described  by  a  Hankel  function  of  the  second  kind  since  the  wave  is  traveling 

radially  outward  and  must  vanish  at  infinity. 

ML  v  m  /  -j  *  &  n  fJ)  /  /)  m  \  -j  h  m 

f  =  f4"  +  4n  c  J  Hn  (#  r)e 

All  wave  functions  will  be  understood  to  have  a  time  dependence  exp, 

Now  consider  a  plane  wave  incident  upon  the  boundary  between  re  cions 
II  and  III.  The  wave  is  described  by  Eq  exp.  -j(K  *  R  -  ot)  where  K  is  the 
propagation  vector  and  R  the  vector  to  a  point  of  observation  (Figure  I). 

Eo  is  the  amplitude  of  the  wave.  The  propagation  vector  may  be  written 
in  the  rectangular  components 


/f'  -  K  A/n  f  £ 

X 

(.») 

f(y  v r.  K s*n  i  £ 

(»*) 

K  =  K  609  i 

(c) 

so  that  E^  exp.  (-j  K  •  R)  is  .1  solution  of  (II). 
Now 


v*  *  K  S  (xto3  5  t  y  S) 

d  —  d  tSf 


»  e 


e#s  (#*f) 


where  r.  ®,  and  z  are  cylincirie.il  co-ordinates.  This  may  be  expanded  into 
cylindrical  wave  functions  (1)  such  that 


.  A  .  «  ,  r  \  -JS#  -J  X*  c*9  ( 

tn  x  EoJ  e' 
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The  wave  function  for  region  III  may  now  be  written  as  the  sum  of  the 
incident  and  the  scattered  wave. 
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VC 


(24) 


From  (21),  (22),  and  (24),  the  separation  constant  h  may  be  inferred  to 
be  the  z  component  of  the  propagation  vector  and  since  this  is  a  rectangular 
component, 


h 


«z 


Kc09  $ 


Then 


j}  =i  iff1—  h*  =  /V  f 

This  solution  for  h  satisfies  boundary  conditions  to  be  imposed  on 
the  fields,  but  at  the  same  time  it  limits  the  aspect  angle  such  that  C  >  O 
for  if  (  is  allowed  to  become  %ero,f^  and  are  not  finite* 

A  more  general  solution  for  h  must  describe  the  propagation  of  energy 
along  the  z  axis  of  the  cylinders*  As  in  the  case  of  circular  waveguides, 
various  modes  of  propagation  are  possible  and  an  infinite  number  of  solu¬ 
tions  of  the  propagation  constant  exist  (4).  If  a  plane  wave  is  incident 
normal  to  the  surface  of  the  cylinder,  no  energy  is  propagated  along  the  z 
axis,  and  (25)  is  an  exact  solution  for  h*  The  remaining  discussion  will 
be  confined  to  this  case* 
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rile  inc  ident  field  may  now  be  considered  to  be  a  superposition 
of  two  lields,  the  firat  having  an  E  vector  parallel  to  the  z  axis  with  its 
associated  H  field  normal  to  z  and  the  second  having  an  E  vector  normal 
to  z  with  H  parallel  to  the  axis.  Then  (|>  is  a  solution  for  E  in  the  first 
case,  and  a  solution  for  H  in  the  second.  The  remaining  components 
may  then  be  found  from  Maxwell*s  equations. 

Since  describes  the  z  component  of  the  fields  and  tangential 
components  must  be  continuous  at  a  boundary,  then  the  wave  functions 
may  be  equated  at  their  respective  boundaries. 

Let  the  incident  field  have  amplitude  Ezo  and  be  polarized  parallel 
to  the  z  axis.  Then 

r*  (A; 

where 


for 

A  =•  K  COS  i  =  O 
From  (1), 

H.  .  -i-  4^- 

9  J€JyU  9r 

Now  equate  the  z  components  of  E  and  the  8  components  of  H  at  the 
boundaries  represented  by  r |  and  r^. 
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From  (<i8)  and  (<J9) 

e„  =  -  A,  i, 

where* 

.;  Wjj)i  [uS'lA.)  1 

Substituting  tor  Cn  in  (26)  and  (27)* 

£„/"*  W„'  (Arz )  -  [W/W-  4,W',)] 

f.. J "  h  OA)]'  ^[^oA.)] 

-  i„  p  [///'(A)]  -  K  i  [TVo]} 
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then 


r f- fs„  Jr [j. (Aj  -_&/ itf *0} 
rf»'  *  B„ '«/*’(/*  O'  3.  £  [AfAJ] 

where 

B„  =  A(A)  -  K  AVAj 


on 


(3^) 


B;  *  a)]'  4»  ^  '*)]  ( 

It  is  more  convenient  to  express  the  derivative  with  respect  to  r  of  the 
Hesse  I  functions  as  the  derivative  with  respect  to  the  argument. 

Thus, 

i  z»  w  -  fi  ip>  *.  w 

Henceforth.  ~7~  -  will  be  written 

d(fir) 


Now 

/ii  (A )  -  (A )  

"  ’■  Bn'  H„<»V>m'z)- A"  V'y  (Aarz) 

where 

B„  =  H ,M (AZ'Z)  -  A*  H*  *  (/ \ ) 

A/r[AYAAr;7A;] 
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The  fields  in  region  III  for  the  parallel  polarization  case  are  now  completely 
defined. 


£  -  K  B„  X  (/>"'*)  *  •».  Hr,  <2>(4mrx  )]  m 

n»-«? 

W  =  "  X  J  *  J”°  [X» *  (/*'*)  +  <Jf» 

n«  -•» 

where  e..  and  e,x  are  unit  vectors. 

4  w 

The  Kirchhoff-Huygens  principle  states  that  if  the  value  of  a  scalai 
field  quantity  is  known  at  every  point  on  any  closed  surface  surrounding  a 
source-free  region,  each  unit  of  surface  can  be  considered  as  a  radiating 
source,  and  the  total  field  at  any  interior  point  is  given  by  integrating  the 
contributions  ot  all  the  individual  elements  over  the  surface.  The  extension 
of  this  principle  to  cover  vector  waves  is  given  by  Stratton  (5)  in  section 
8.  14.  We  consider  the  field  at  an  interior  point  (P)  of  a  volume  bounded 
by  the  cylinder  and  infinity. 

The  electric  field  at  P  is  given  by 

E=  J  j  -(n-i)v+  d*,  H0) 

3 

where  n  in  a  unit  vector  normal  to  the  surface  and  the  field  quantities  in 


I  t 


the  integrand  are  those  just  inside  the  surface  enclosing  the  volume.  Since 
the  fields  vanish  at  infinity,  the  quantities  in  the  integrand  are  defined  by 
equations  (38)  and  (39)  if  £  is  replaced  by  £^.  Phi  is  the  Green's  function 
for  free  space  defined  by: 


where  R'  is  the  distance  from  a  point  on  the  surface  to  the  point  P. 

For  large  R',  relative  to  the  wake  length  the  gradient  may  be 


written 


and  the  1/R'  attenuation  factor  may  be  approximated  by  1/R0,  where 
R0  is  the  normal  distance  from  the  cylinder  axis  to  the  point  P.  However, 
a  more  exact  value  is  necessary  for  the  phase  factor  and  this  will  be 

R  »  +  Yz  COi  6 


The  Rq  in  the  phase  factor  may  be  neglected  since  it  is  constant.  The 


Green's  function  and  its  gradient  are  now  defined  by 

J  K  C 


The  geometry  is  illustrated  in  Figure  2.  Replacing  the  vector  quantities 
by  the  appropriate  scalars,  equation  (40)  may  now  be  written: 
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where  it»  the  length  ui  the  cylinder. 


r2K'i 
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(41) 


Tilt-  integral:*  in  (41)  >irv  tin*  integral  representations  of  IWssel  iunctmns. 


Z  tr  j'n  Jn(K'*)  «  f  *  ■*'•***'**"•  da 
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Substituting  (4£)  and  (41).  (41)  becomes 


IS 


Figure  l 

>os»-Sertion  Geometry 
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J'he  function  in  brackets  is  the  Wronskiun  of  the  Bessel  functions  (2)  and 


is  given  by 


ill 

Jn  (/  r. 


)  w 


H„  (/  / 


Therefore, 


^  t4r0£**Yi.(~ifi*n  (44> 

When  the  incident  field  is  polarized  transverse  to  the  z  axis,  the 
same  method  is  used  to  calculate  the  magnetic  field  which  is  given  by: 
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TTK0 
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where 

J*  7/4mr%)  -  £„ '  J, 


K  -  n.0,(/t\y  o„ 


(45) 
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(4H) 
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(49) 
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A‘H„W($‘r,)  J.‘ (A.) -/*  K™' (/>*',)& f/1',) 


D„  = 


Tin*  ii  field  is  (lion  found  .it  P  by 
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H  -=  “  (  | :J6>t(n*£)+  ■■  (A  *H)*V+  (h-H)Vt j>]  da 


where  till  quantit ies  are  as  previously  defined* 

For  arbitrary  polarization,  the  amplitudes  ot  the  incident  fields  h 
and  llzo  arc  related  by  (he  polarization  angle. 

The  radar  cross-section  crr  is  defined  by 


4  ir  H 


fT 

jv1 


(SO) 


where  the  N's  are  time  averaged,  received  and  transmitted  Foynting  vectors 
and  R  is  the*  distance  to  the  target,  and  equal  to  K^,  Since  N  is  proportional 
to  the  square  of  the  fields,  we  may  write 


ar  =  MrrA V 


£f- 


4 rr 


*|  H 


*xZ 


For  parallel  polarization  (SO)  becomes 

1  I  WF*  '  •  -  *  2 


£t 

tr 


!H 


For  traiiavt-rHc  pul.tr i/ation  (50)  becomes 

Z  ...  ..  .2, 
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The  solution  to  equations  (SI)  and  (S^i)  has  been  pro^ramimd  tor  the 
IBM  7090  computer.  The  initial  input  data  are  determined  .ns  accurately 
as  possible  front  flow  field  analysis  (7).  These  data  are  then  automatically 
incremented  to  cover  all  values  assumed  to  exist  in  a  laminar  wake.  l*he 
tallowing  method  is  used  to  obtain  the  data. 

1.  An  arbitrary  value,  near  unity,  is  assumed  for  the  ratio*  coju)^ 
in  the  inner  cylinder. 

1.  The  largest  radius  for  which  condition  one  holds  is  determined 
from  flow  field  calculations  and  becomes  r^. 

i.  The  radius  is  determined  from  How  field  calculation- 
such  that  is  equal  to  two  at  the  outer  boundary. 

4.  The  complex  permittivities  for  regions  one  and  two  arc  a\ cram's 
of  the  values  found  between  the  boundaries. 

5.  The  complex  arguments  of  the  Bessel  functions  are  incremented 
systematically  by  me  remolding  the  modulus  and  the  argument. 

An  attempt  is  being  made  to  keep  the  increments  small  enough  to 
observe  resonances  and  yet  vary  the  parameters  sufficiently  to  co\cr  all 
possible  situations.  Results  of  these  computations  v/ill  be  reported  at  a 
later  date. 
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